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Abstract 

Monte Carlo simulations based on a first-principles-derived Hamiltonian are conducted to study 
the properties of Pb(Zri_ a; Ti ; j;)03 (PZT) alloys compositionally modulated along the [100] pseu- 
docubic direction near the morphotropic phase boundary (MPB). It is shown that compositional 
modulation causes the polarization to continuously rotate away from the modulation direction, 
resulting in the unexpected triclinic and C-type monoclinic ground states and huge enhancement 
of electromechanical responses (the peak of piezoelectric coefficient is as high as 30000 pC/N). 
The orientation dependence of dipole-dipole interaction in modulated structure is revealed as the 
microscopic mechanism to be responsible for these anomalies. 

PACS numbers: 77.80.-e, 77.84.Dy, 77.22. Gm 
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Complex perovskite alloys have been studied for more than fifty years, which are of great 
importance not only because they have wildly technical applications for their large elec- 
tromechanical responses, [l] but also because they are fundamentally interesting for their 
anomalous properties and the relations to the microscopic structure. Recently, two break- 
throughs were achieved in the research of complex perovskites. One is the observation of large 
enhancement of piezoelectric constants in single crystals. The second is the discovery of 
a low-temperature monoclinic phase in Pb(Zri_ x Ti x )03 (PZT) near the morphotropic phase 
boundary (MPB).0, Q] Inspired by such progresses, lots of works were done to search larger 
electromechanical response and unobserved phases in complex perovskites. Particularly, an 
ab initio calculation on Pb(Sco.5Nbo.s)03 (PSN) revealed that atomic rearrangement would 
lead simultaneously to large electromechanical responses and unusual structural phases. Q] 
The difference between the valences of the B atoms (Sc 3+ and Nb 5+ ) was pointed out to 
be the main reason for the existence of these anomalous properties and this mechanism is 
expected to be generally applicable in a given class of heterovalent perovskite alloys, which 
is likely to have large technological and fundamental implications. However, it is unclear 
whether there is any interesting effect in the modulated homovalent systems where the B 
atoms have the same valences. 

In this letter, we investigate the influences of compositional modulation on the homovalent 
PZT solid solutions near MPB. Remarkably, it is found that the atomic ordering along the 
[100] direction leads to the suppression of P x (x-component of the polarization). Unexpected 
new ground states of triclinic (Tri) and C-type monoclinic (M<?) phases appear with greatly 
enhanced piezoelectric coefficients. The orientation dependence of dipole-dipole interaction 
in modulated structure is revealed to be responsible for these anomalies. Such microscopic 
mechanism is expected to work in MPB region of any homovalent perovskite when the atomic 
ordering is produced. 

We adopt the effective Hamiltonian of PZT alloys proposed by Bellaiche, Garcia and 
Vanderbilt, 9] which is derived from first-principles calculation, to predict the properties of 
compositionally modulated PZT structures by conducting Monte Carlo simulations. In this 
scheme, the total energy E is written as the sum of an average energy and a local energy 
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+E loc ({u i },{v i },{a j }), (1) 

where Uj is the local soft mode in unit cell i, v, is the dimensionless local displacement 
related to the inhomogeneous strain^], i]h is the homogeneous strain tensors and o~j = 
±1 represents the presence of a Zr or Ti atom, respectively, at lattice site j of the PZT 
alloy. All the parameters of Eq. (1) are derived from the first principle calculations and are 
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listed in references. |9J, HO] The modulated structures under consideration are made of the 
sequences Pb(Zri_ a;+l ,Ti x _ ! ,)03/Pb(Zri_; E _j,Ti. r+l ,)03 along the [100] direction with B atoms 
randomly distributing within each [100] plane. 10x10x10 or 12x12x12 supercells with 
periodic boundary conditions are used in Monte Carlo simulations to get well converged 
results. Eq. (1) is based on the virtual crystal approximation (VCA).^, Q] The parameters 
of Zr and Ti atoms are composition-dependent, which reflects an average influence from 
surrounding atoms. Considering that the modulation length of our studied structures is 
small (only two monolayers) and a Zr or Ti atom will be influenced by atoms from planes with 
different compositions, we adopt the parameters corresponding to the overall composition x 
for all atoms in the modulated structures. Some detailed tests are also carried out, which 
show that the anomalous properties of the modulated structures revealed in this paper are 
not sensitive to the choice of parameters. 

First we consider the properties of modulated PZT as a function of composition x by set- 
ting parameter v to its maximum value {y — x). Fig. 1(a) shows the cartesian coordinates 
(u x , u y , u z ) of the supercell average of the local mode vectors at T=50K. (Due to the differ- 
ence between the theoretical and experimental temperatures,^ it corresponds to an "exper- 
imental" temperature around 30K.) In the non-modulated case, it was demonstrated that 
disordered Pb(Zr 1 _ a; Ti 3 .)03 has a pseudocubic [111] direction polarized (u x = u y = u z > 0) 
rhombohedral phase (R) for x < 0.475, a [001] direction polarized (u x = u y = < u z ) tetrag- 
onal phase (T) for x > 0.49, and a A-type monoclinic phase (M^) with < u x = u y < u z 
for 0.475 < x < 0.49. In the modulated case, for compositions outside MPB, the local 
mode vectors are slightly affected as shown in Fig. 1(a). However, the precise symmetry has 
been changed by the atomic ordering, and the phases are recognized as B-type monoclinic 
(Mb) and orthorhombic (O). The interesting phenomena appear within the MPB region. 
For 0.480 < x < 0.485, u x decreases drastically, which results in a triclinic (Tri) phase with 
Ux 7^ Uy 7^ u z . When x is larger than 48.5%, u x vanishes and an unexpected ground state 
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FIG. 1: (a) Average cartesian coordinates of the local mode vectors; (b) piezoelectric coefficients 
di5 and g?24! (c) dielectric susceptibilities xn anci X22 as functions of composition x in atomic- 
ordered PZT alloys made of the Pb(Zri_ a;+i ,TL E _ i ,)03/Pb(Zri_ x _ I ,Ti :c _|_ I ,)C)3 sequences, v is set to 
its maximal value as v = x. The temperature in simulation is 50K. 

appears in the modulated systems: a C-type monoclinic (Mc) phase with = u x < u y < u z . 
Fig. 1(a) thus depicts a polarization rotation in modulated PZT structures following the 
sequence of M^-Tri-Mc-O as the Ti composition increases. Although Tri and Mc phases 
in PZT were first demonstrated in the simulation of electric-field induced transformation in 
disordered rhombohedral PZT by Bellaiche et a/.,^^ it has not been observed experimen- 
tally yet because the predicted electric field required to induce the Tri and Mc phases is too 
high (~ 500 kV/cm for composition close to MPB) compared with what experiments can 
provide. Thus the present study demonstrates an alternative routine to observe such phases 
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FIG. 2: (a) Average cartesian coordinates of the local mode, and (b) piezoelectric coefficient di§, 
as functions of the parameter v in modulated PZT structures with x = 0.485 at 50K. 

with less experimental obstacles. 

Fig. 1(b) shows the corresponding piezoelectric coefficients and Fig. 1(c ) the dielectric 
susceptibilities which are calculated in the correlation-function approach [13(. It shows that 
d\5 reaches its maximum near the Tri to Mc transition and remains a remarkably large value 
above 1000 pC/N in MPB region, while ^4 achieves values above lOOOpC/N in Mc phase 
and reaches its maximum near the Mc to O transition. Another important piezoelectric 
coefficient, ^33, is not shown here because the modulations scarcely affect u z and have little 
effect on <i 33 . The existences of the M^-Tri, Tri-Mc and M^-0 transitions also result in the 
peaks of dielectric responses, as shown in Fig. 1(c). 

We also investigate the effects of the modulation parameter u, the fluctuation amplitude, 
on the structural and piezoelectric properties of the modulated structures which have overall 
composition x = 0.485. Fig. 2(a) shows the average local modes as a function of v at 50K. 
The structure of v = (in the case of disordered PZT alloy) has u x = u y < u z character- 
izing the Ma structure with polarization lying between [111] and [001] directions. When v 
becomes larger than 0.155, u x begins to decrease. With increasing u, u x decreases gradually 
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FIG. 3: Piezoelectric coefficients in modulated PZT structure with x = 0.485 and v = 0.305 as 
functions of temperature. 

whereas u y and u z increase very slowly. This behavior corresponds to the phase with a Tri 
symmetry. When v becomes larger than 0.45, u x becomes null, while u y and u z remain 
unequal, indicating a Mc symmetry with a polarization between [011] and [001] directions. 
Summarily, increasing the parameter v leads to a continuous rotation of the polarization 
and transformation between different ground states which are driven by the annihilation of 
u x — the polarization component along the direction of compositional modulation. 

Theoretical and experimental studies have suggested that polarization rotation is essential 
for the large piezoelectric response in the perovskite alloys. Q| Such an effect is clearly 
demonstrated in Fig. 2(b) where the piezoelectric coefficient di§ is depicted. With decreasing 
u x , di^ increases from 2000 pC/N to 8000 pC/N gradually in the Tri phase for 0.2 < v < 0.4. 
Then dis reaches a very steep peak about 30000 pC/N at v ~ 0.4 where the polarization 
most easily rotates [see also Fig. 2(a)]. After that, d\^ decreases slowly but still remains 
large values even in M c phase. As one can see in Fig. 2(b), d\*> is larger than 10000 pC/N 
for a broad range of v > 0.4. To our best knowledge, such a huge enhancement of the 
piezoelectric coefficient has never been reported in any perovskite material before. 

Temperature properties of the material are important for its application performances. 
Here we investigate the influence of temperature on a structure with x = 0.485 and v = 0.305, 
which adopts a Tri ground state in Fig. 2(a). The resulting piezoelectric properties are 
depicted in Fig. 3. A huge d 2 4 piezoelectric coefficient peak appears near the O to Mq 
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transition (the information of average local mode is not shown here), around which c?24 
remains larger than 5000 pC/N for a broad range of temperature. The simulations also 
predict the d\*, peak consistent with the transition from Mq to Tri phase. In addition, it is 
found that with increasing u, the M^-Tri transition temperature decreases, and consequently, 
the temperature of c? 15 peak is dependent on the value of the modulation parameter v. This 
dependence brings great convenience to developing piezoelectric devices for application at 
different temperatures. 

The above intriguing results can be understood by considering the orientation dependence 
of dipole-dipole interaction in modulated structures. 

Although PbZrC>3 and PbTiC>3 are similar enough to form Pb(Zr 1 _ x Ti I .)03 solid solutions 
at any composition x, their intrinsic polarization properties are different to some extent. As 
a result, the magnitude of polarization vibrates with varying local composition that is caused 
by the structural modulation. In our calculation, the Ti-rich planes have larger average local 
mode, while Ti-poor planes have smaller values of local mode. 

In the homovalent alloys such as PZT, the long-range polarization interactions are de- 
scribed in dipole-dipole interaction energy: U| 



_Z ^ Uj • Uf - 3(Ry • Ui)(Rij • Uj) . . 

^dipole - — 2^ ' K Z ) 

too {<j H-ij 

where Z* is the Born effective charge and the optical dielectric constant of the material. 
When the dipole is uniform = (u)), -Edipoie is isotropic, i.e., .Edipoie is independent of 
the direction of the polarization. However, when the magnitude of the dipole fluctuates due 
to the structural modulation, the situation differs. -Edipoie becomes dependent on the angle 
between the directions of modulation and polarization [see Fig. 4(a)]: when the polarization 
is parallel to the modulation, the dipoles are aligned shoulder to shoulder and the energy 
is higher; when the polarization is perpendicular to the modulation, the dipoles are aligned 
end to end so that the energy is lower. To further clarify this point, we present a simple 
numerical analysis here. The local modes in a [100] plane are assumed to fluctuate as 

u i = (u)[l + ZK i ), (3) 

where £ is the fluctuation amplitude, and Ki equals to +1 or —1 depending on which [100] 
plane the dipole locates. It is easy to show that the dipole-dipole interaction energy can be 
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FIG. 4: (a) Schematic graphics of dipole alignment in modulated structure: dipole aligned parallel 
to modulation direction with higher energy (the left), and alignment perpendicular to modulation 
with lower energy (the right), (b) Variation of dipole-dipole interaction, AE^ipoic [see Eq. (4)], 
as a function of ip (the angle between the polarization and the modulation direction). AE^ipoie * s 
measured in units of -Ey ipole , while tp in units of tt/2. 



rewritten as 



-^dipole — -^dipole + £ 2 A-Erfipole- 



(4) 



A-Edipoie reflects the variation of dipole energy caused by the structural modulation and re- 
lies on the angle between the polarization and the modulation directions (denoted as ip). In 
Fig. 4(b), we plot A£y ipo i e as a function of tp with polarization rotating from [100] to [010] 
direction. It shows that AE^ipoie nas its maximum when the polarization is parallel to the 
modulation direction. When the direction of polarization deviates from that of modulation, 
AEdipoic continuously decreases and reaches the minimum with polarization perpendicular 
to the modulation direction. Accordingly, compositional modulation will cause the polar- 
ization to rotate away from the modulation direction, which results in unusual structure 
transitions and induces huge enhancement of piezoelectric coefficient, as revealed in the 
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above case of PZT. In perovskite ultrathin films, the polarization is restricted in the normal 
direction by the strain from matrix, so the dipole-dipole interaction increases and violates 



the ferroelectricity below a critical thickness. 



16] 



It is noted that the structures outside MPB region are not sensitive to the atomic ordering. 
The main reason is that the dipole-dipole interaction is not the only part to determine the 
Hamiltonian. When the rotation occurs, the other parts of the Hamiltonian may increase 
and cause an energy barrier. MPB is a region where different phases have close free energy, 
which means that the energy barrier of polarization rotation is small. So the dipole-dipole 
energy dominates in this region and produces the anomalies. When the composition deviates 
from MPB, such energy barrier increases and the polarization rotation is prevented. 

On the basis of the general microscopic mechanism discussed here, it is expected that 
solid solution made of homovalent perovskite ferroelectrics with a MPB should have the 
similar structural and piezoelectric anomalies as those of PZT when the atomic ordering is 
arranged. Besides, compared with heterovalent alloys such as PSN, experimental modifica- 
tions to the atomic ordering of homovalent alloys such as PZT is more practicable because 
Pb(Zri_ a Ti a; )03 is stoichimetric at any x value. 
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